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Abstract. - Energy dissipation is studied for a hard magnetic tip that scans a soft magnetic 
substrate. The dynamics of the atomic moments are simulated by solving the Landau-Lifshitz- 
Gilbert (LLG) equation numerically. The local energy currents are analysed for the case of a 
Heisenberg spin chain taken as substrate. This leads to an explanation for the velocity dependence 
of the friction force: The non-linear contribution for high velocities can be attributed to a spin 
wave front pushed by the tip along the substrate. 



I 

o 
o 



(N 
> 

O 

o 



> 



Introduction. — On the macroscopic scale the phe- 
nomenology of friction is well-known. However, investi- 
gations of energy dissipation on the micron and nanome- 
ter scale have led in recent years to many new insights 
[l]. This progress was made possible by the development 
of modern surface science methods, in particular Atomic 
brce Microscopy, which allows to measure energy dissi- 
pation caused by relative motion of a tip with respect to 
a substrate. 

Studies concerning the contribution of magnetic degrees 
of freedom to energy dissipation f^J^ form a young sub- 
field of nanotribology, which has been attracting increas- 
ing interest in recent years. Two classes of models have 
been considered, which show different phenomena. The 
first one is Ising-like spin systems with two equivalent half 
spaces moving relative to each other [4j|8]. In this case, 
friction is induced by thermal fluctuations, and hence is 
not present at zero temperature. In the second class of 
models [9 - 13 , there is no symmetry between slider and 
substrate: The slider, representing e.g. the tip of a Mag- 
netic Force Microscope, interacts only locally with a pla- 
nar magnetic surface. While scanning the surface, the tip 
in general excites substrate spins and hence experiences 
friction, even at zero temperature. The present study be- 
longs to the second class of models. 

We investigate the nature of the substrate excitations 
caused by the tip motion for a classical Heisenberg model 
with Landau-Lifshitz-Gilbert (LLG, [MjjTs]) dynamics 
(precession around, and relaxation into the local field di- 



rection). As the spins are continuous variables, spin wave 
excitations are possible. As we will show in the following, 
their properties are reflected in the velocity dependence of 
the friction force. Spin waves are increasingly attracting 
interest: e.g. in the last years a new subfleld of magnetism, 
magnonics, has been developed, where materials are stud- 
ied with respect to their spin wave properties 16p7 . One 
motivation is to create new devices using spin wave logics 
or novel concepts of data storage. 

In a previous work we showed that friction in this model 
is proportional to the scanning velocity v ("viscous be- 
haviour"), provided that the tip does not move too fast 
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10 . The reason can be found in continuous excitations, 
while the motion in the Ising-model consists of discrete ex- 
citations and relaxations, which yields a constant friction 
force for low v. In the present paper we focus on the lo- 
cal dissipation processes in order to explain, why for high 
velocities deviations from the viscous behaviour exist. 

Simulation model. — To simulate a solid magnetic 
material, we consider a chain of N~320 classical, nor- 
malised dipole moments ("spins", cf. fig. [T]) Si—fj,^/fj,s, 
where /i^ denotes the material-dependent atomic magnetic 
moment. The spins represent magnetic moments of single 
atoms, arranged with a lattice constant a along the x-axis. 
Two lattice constants above the spin chain, a magnetic tip 
Stip moves with constant velocity v—ve^. Its magnetisa- 
tion is fixed in z— direction. At the beginning of each 
simulation the tip is positioned at the centre of the chain. 

In order to keep boundary effects small, we use a con- 
veyer belt technique with anti-periodic boundary condi- 
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Fig. 1; Snapshot of a simulation. The colour encoding denotes the spins' orientation in the i/z-plane (the tip is moved along 
the spin chain axis to the right). In front of the tip a spin wave is visible, i.e. an oscillation of the spins around the x-axis. 



tions: When the tip has moved by one lattice constant, 
the boundary spin at the back end is deleted and a new 
spin with opposite direction is added at the front end of 
the chain. This shift puts the tip back to the centre of the 
simulation cell. 

This paper analyses, how local excitations contribute 
to magnetic friction. Each substrate spin contributes its 
exchange interaction, Hg^i^, and its interaction with the 
tip field, Hjip to the Hamiltonian 
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As the substrate spins represent a ferromagnetic solid, we 
use the anisotropic Heisenberg-Hamiltonian, 



n. 



(i) 
sub 



J 



Si • (Si+1 + Si_i) — dzSf 



(2) 



J>0 describes the ferromagnetic exchange interaction be- 
tween Si and its nearest neighbours i±l. In order to avoid 
double counting, half of the pair interaction is attributed 
to either spin. dz——0.1J is the anisotropy constant: Here 
it defines an easy plane anisotropy, thus the substrate 
spins prefer an alignment in the xy-plane. The moving 
tip interacts with each substrate spin by the dipolar inter- 
action. 



(i) _ _ 3(Si ■ ei)(Stip ■ e.t) — S,; ■ Stip 
tip - ' 



where Ri — |Ri| is the length of the distance vector = 
Tj — Ttip, and Bi its unit vector = Hi/Ri. and rtip 
denote the position vectors of the substrate spins and the 
tip, respectively, w quantifies the dipole-dipole coupling 
of the substrate and the tip, with w |Stip| = 10 Ja^ in this 
paper. 

While the tip magnetisation direction is fixed in time, 
the substrate spins are allowed to change their orientation. 
To simulate their dynamics, we solve the LLG equation. 



7[S., X h, + aS^ X (S, x h^)] 



(4) 



numerically via the Heun integration scheme, where 7 = 
7 [/Lts(l -|- a^)]~^ with the gyromagnetic ratio 7. The first 
term represents the Lamor precession of each spin in the 
effective field, 



with the precession frequency 7|hi|. The precessional mo- 
tion preserves energy. Dissipation is introduced by the 
second term which causes an alignment towards h, . a is a 
material constant which can be obtained from ferromag- 
netic resonance experiments and represents the coupling of 
each spin to a reservoir of zero temperature. By adding a 
stochastic term to the effective field, it is possible to study 
the influence of finite temperatures as done in 
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However, in order to analyse the non-equilibrium excita- 
tions it is advantageous to suppress thermal spin waves by 
setting temperature equal to zero in this work. 

In order to discuss frictional losses occurring in the sys- 
tem the global energy balance was analysed in 10 : 



~ pump J diss 7 

P -VpW - V "tip 

,^1 "^tip 
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^i--E^iis-E^"(s.xh,r 



(6) 
(7) 

(8) 



The only explicit time-dependence of the Hamiltonian H 
stems from the motion of the tip. It leads to the first term 
in eq. (|6| , which is the energy pumped into the spin system 
per unit time by an outside energy source that keeps the 
tip moving. Accordingly we call it the "pumping power" . 
(3) Its local contribution, Ppumpi is the energy transferred per 
unit time from the tip to substrate spin S^. The friction 
force, F = — Fe-c, the substrate exerts on the tip is given 

by 



F 



(Ppump) 



(9) 



where the angular brackets denote a time average over at 
least one period a/v. 

^diss represents the energy current from spin Sj into the 
heat bath. In other words, this is the energy dissipated at 
site i per unit time. Dissipation always occurs when the 
system relaxes towards the ground-state, in which the spin 
at site i is aligned with the local field-direction h^. With- 
out tip movement, Ppump is zero and Pdiss leads the system 
quickly to its ground state. For a tip moving at constant 
velocity, a steady non-equilibrium state is reached, where 
the time averaged derivative {dH/dt) vanishes, because 
all energy pumped into the system is dissipated. Then 
the two power terms in cq. ^ cancel. 
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Fig. 2: Illustration of the local energy balance. The arrows rep- 
resent the directions in which the energy transfers are counted 
positive: a positive j'^' rises the energy at site (i), but lowers 
the one at site (i— 1). 



When evaluating the local energy balance instead of 
eq. ([g]), energy currents je within the substrate have to be 
taken into account, which transport energy from one spin 
to its neighbour (cf. fig. [2]). By taking the time derivative 
of the local Hamiltonian one obtains: 



^^^^^ =J^lp-a(divjE)(')-/^:i> 



dt 



- - J(S, - S,_i) 



diss ' 
Si-1 + Si 



J 



(Si • Si_i — Si_i • Si) 



(10) 
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Simulation Results. — Let us consider the steady 
state in a co-moving frame: the local quantities do not 
depend on spin index i and time t separately, but only 
on the (continuous) coordinate Xi = Ri(t) • e^. For con- 
siderations, where all spins are equivalent, we can drop 
the index i, e.g. the tip position is always at a; = 0. 
In its vicinity, fig. [sjja) shows the local pumping power, 
as well as the discrete divergence of the energy current, 

(divjE)^*^ = ~jE^) /'^^ ^ functions of x. 

The physical interpretation of fig. [sj^a) is the following: 
When the tip approaches, a substrate spin lowers its en- 
ergy by adjusting to the inhomogeneous tip field at the 
cost of the exchange interaction. When the tip has passed 
by, it returns asymptotically to its higher energy in the 
absence of the tip field. This means that the tip injects 
energy Ppump(2;) oc v per unit time at x < and extracts 
apparently the same amount of energy from the substrate 
spins at X > 0. With respect to origin and curve shape, 
this is very similar to an electrical charge passing by a 
charge of opposite sign on a straight line. The apparent 
central symmetry holds only up to first order in v, though. 
The small asymmetry, not noticeable in figjsja), is due to 
dissipation, which will be discussed below. But first we 
derive the steady state current within the chain. 

In the steady state, we have 



Si = S(a;i) = vdxS{xt) = v 



Si+i — Si 



(12) 
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Fig. 3: (a) Pumping power and transported power, rescaled 
by velocity, at a = 0.1. Negative pumping power represents an 
energy return from the spin chain to the tip, positive an energy 
injection from the tip into the chain, (b) Energy current for 
five different time steps (points, for a = 0.1), as well es the 
dissipation power (solid lines). They turn out to be exactly 
proportional to each other, the coefficient being linear in av. 



tient, holds up to first order in a. Plugging that into eq. 



the current jg"* reads 

7«-— fl- 



Si-i • Si+i) 



(13) 



This shows that je transports the exchange energy to be 
paid for orienting the spins according to the inhomoge- 
neous tip field. Correspondingly, the source of this current 
(divjE > 0) is behind, and its sink (divjE < 0) is in front 
of the tip, as seen in fig. [sjja). 

The terms discussed so far are reversible and hence in- 
dependent of the damping constant a: To first order in v 
they add up to zero in eq. ([6]). The origin of dissipation 
is that the spin pattern does not follow the tip instanta- 
neously, but with a delay, which corresponds in the co- 
moving frame to a lag Ax on av 10 . It is a manifestation 



of the driving out of equilibrium, which the spin relaxation 
must counteract. Fig. |3jb) shows that Pdiss(a;)/jE(a;) is 
indeed a constant cx av, which we may call the driving 
forc^ As pointed out in fig. [sjja), je is proportional to 
V. This only holds true for velocities not much larger than 
Vq « 0. 3I7 Ja//is (see eq. ( [T6| below). This imphes that 
for small velocities Pdiss c>c mr , which gives rise to a fric- 
tion force, eq. proportional to av. 



where the third equality, due to using the difference quo- 



^This is analogous to electrical power P = UI, where the voltage 
U provides the driving force for the current I. 
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Fig. |3jb) shows an important qualitative difference be- 
tween energy currents for tip velocities below, respectively 
above vq. For low velocities, the energy current is concen- 
trated around the tip position in an essentially symmetric 
way: Whatever exchange energy is released behind the tip, 
is reabsorbed in front of it. For high velocities, however, 
an additional shoulder in front of the tip appears. This 
shoulder represents a part of the energy current, which 
can leave the tip's immediate neighbourhood and propa- 
gates further along the spin chain, until it is damped out. 
We call this contribution non- confined. The propagation 
range depends on the damping constant a, as can be seen 
in fig.|4] The lower the damping constant, the farther the 
current extends. 

In order to evaluate this quantitatively, we define the 
non-confined energy current as 



jnc{x) ^ i^{x) - j^{-x) for x>0. 



(14) 



It is plotted for several a- values in fig. |4]jb). The axes 
are rescaled in order to show that the range shrinks 
with increasing damping approximately like q;~°'^, and 
that the amplitude of the non-confined current also de- 
creases roughly like a~°'*. The range and the ampli- 
tude of the non-confined current combine in such a way, 
that the integral over jnc{x) is nearly proportional to a~^. 
Hence, when multiplied by the driving force oc au, the a- 
dependence nearly cancels. The contribution of the non- 
confined excitations to friction is therefore approximately 
independent of a, in contrast to the confined contribu- 
tion discussed above. As will be explained below, the two 
contributions also depend differently on velocity. 

The non-confined excitations can be regarded as spin 
waves. An excitation means a deflection of a spin Si from 
the local field , leading to dissipation at the correspond- 
ing site according to eq. ([8| and a precession around in 
a plane perpendicular to h^. Accordingly, 



hi X Gy 

|h,:xe„ 



and e,;,f, = 



hiXe,^a 
Ih-jXe^.a! 



(15) 



form an appropriate local basis to illustrate the excita- 
tions. Spins far in front of the tip experience a field which 
points in x-direction, thus the ei^a-component points in 
z-direction. Near the tip the basis changes as sketched in 
the inset of fig. [Sj 

In the low velocity regime, a deflection from the local 
fields is present solely in the vicinity of the tip, according 
to the purely confined contribution to friction discussed 
above. For velocities, where non-confined currents can 
be observed, additional oscillations are present in front of 
the tip. A Fourier analysis of our simulation data shows 
that their wavenumber has approximately a linear velocity 
dependence: 



oc w— Wo wi 



th 



0.31 ^Jaj [Is 



(16) 



The resulting empirical coefficient 0.31 must be expected 
to depend on system parameters like the tip field's shape 
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Fig. 4: (a) Energy currents for «=0.6 and some a values. When 
the damping is lowered, the current can proceed further in the 
substrate, (b) Rescaled non-confined current vs. a rescaled 
position for w=0.85. 



and amplitude, which determine the spin waves' confine- 
ment. For small enough damping, the oscillations presum- 
ably extend arbitrarily far in front of the tip. Asymptoti- 
cally we can neglect the tip field and consider an isolated 
spin chain with exchange interaction only. In its ground 
state, all spins point into the S8-nie direction, Sciy Gx- 

For 

small perturbations of S; from this direction, the LLG- 
equation with 

h,=J(S,_i+S,+i) (17) 

can be linearised (e.g. [l8]). The solution is a spin wave 
with an oscillating part of 

Si — By S cos{ika — cot) + Bz S sin(zfca — cot) (18) 

to first order in its small amplitude d. Its dispersion rela- 
tion 



u}(k) = 4 — ^(1- cos(fca)) 

Ms 



(19) 



yields a group velocity w oc fc in the long wavelength limit. 



The finding eq. ( 16 ) indicates that this holds true even in 
the more complicated system with the inhomogeneous tip 
field. 

Inserting Si — + Si into eqs. (17) and (Isl) yields for 



the wave of wave number k a dissipation of 

A I ka^ 

Pdiss(fc) oc sin 



(20) 



Using eq. ( 16 1 and assuming that the amplitude of the spin 



wave excitations 5 does not depend on v and that fca^l. 
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Fig. 5: Spin component perpendicular to the corresponding 
field. While for low velocities one precession around the local 
field is observable, for higher velocities more precession cycles 
in front of the tip are present. For lower damping excitations 
reach further. Sketch: Definition of the field dependent basis. 
For two sites spin and field values are sketched. For 5*2 the 
local field points in x— direction (as the situation in the studied 
system in front of the tip is) , and the spin precesses on the blue 
disk. The shift of the field for Si to the bottom (as it may be 
induced by the tip-field) yields also a change of the disk and 
the appropriate basis. 



this predicts a non-linear velocity dependence of the spin 
wave contribution to friction like ^"-^ . 

V 

The total magnetic friction force is thus predicted to be 
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Fig. 6: (a) Friction force for several velocities, with the non- 
confined part marked by the shading, (b) Non-confined con- 
tributions to the friction force times velocity as well as the 
function B(v — uo)*, where vq has been taken from eq. ( 16 1. 



open question. Studies dealing with this are already in 
progress and will be reported in a future work. 



F ^ Aav + B{a)Q(y - vq)- 



(21) 



where Q(x) is the Heaviside step function, and the coeffi- 
cients B as well as vq may depend on system parameters 
like the tip field. In fig. [6] this total force is plotted, and 
the simulation results are in good agreement with eq. (21 ). 



Conclusion and outlook. In this work, wc could 
separate two distinct contributions to magnetic friction by 
examining the energy current in a spin chain. The con- 
fined current results in a friction force F=Aav, which is 
in accord with our earlier results for 2d [lO] and 3d [l2] 
substrates. Above a threshold velocity vq, spin wave exci- 
tations may leave the tip's immediate neighbourhood and 
form a damped wave packet in front of the tip, propagat- 
ing along with it. These excitations are the stronger, the 
weaker the damping. They lead to an additional contri- 
bution to friction with a non-linear velocity dependence. 
The dependence of the non-confined contribution on a is 
not trivial, because the range as well as the amplitude of 
the energy current are influenced in a non- linear way, cf. 
figgb). 

Important extensions of the present investigation in- 
clude the influence of dimensionality on the non-confined 
spin waves. Here the propagation is not confined to the 
tip's motion direction. The influence of thermal spin waves 
and their interaction with the free spin waves is another 
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